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We present inclusive charged hadron elliptic flow (^2) measured over the pseudorapidity range 
\ri\ < 0.35 in Au+Au collisions at ys„„ = 200 GeV. Results for V2 are presented over a broad range 
of transverse momentum {pr = 0.2-8.0 GeV/c) and centrality (0-60%). In order to study non-flow 
effects that are not correlated with the reaction plane, as well as the fluctuations of V2 , we compare 
two different analysis methods: (1) event plane method from two independent sub-detectors at 
forward (I77I = 3.1-3.9) and beam (I?;! > 6.5) pseudorapidities and (2) two-particle cumulant method 
extracted using correlations between particles detected at midrapidity. The two event-plane results 
are consistent within systematic uncertainties over the measured pr and in centrality 0-40%. There 
is at most 20% difference of the 112 between the two event plane methods in peripheral (40-60%) 
collisions. The comparisons between the two-particle cumulant results and the standard event plane 
measurements are discussed. 
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PACS numbers: 25.75.Ld 



I. INTRODUCTION 

Collisions of Au+Au nuclei at the Relativistic Heavy 
Ion Collider (RHIC) produce matter at very high en- 
ergy density [11, d, d, [!]• The dynamical evolution of 
this hot and dense medium reflects its state and the de- 
grees of freedom that govern the different stages it under- 
goes 0, 0, 0] ■ Azimuthal anisotropy measurements serve 
as a probe of the degree of thermalization, transport co- 
efficients and the equation of state (EOS) [81,1^,1101 of the 
produced medium. 

Azimuthal correlation measurements in Au-|-Au colli- 
sions at RHIC have been shown to consist of a mixture of 
jet and harmonic contributions [ll|, [H, [3| • Jet con- 
tributions are found to be relatively small for pr ^ 2.0 
GcV/c, with away-side jet yields strongly suppressed p^ . 
Significant modifications to the away-side jet topology 
have also been reported [3, [H, [O] ■ The harmonic con- 
tributions are typically characterized by the Fourier co- 
efficients, 

w„ = (cos(n[0-«'Rp])) (n = l,2,...), (1) 

where represents the azimuthal emission angle of a 
charged hadron and is the azimuth of the reaction 
plane defined as containing both the direction of the im- 
pact parameter vector and the beam axis. The brack- 
ets denote statistical averaging over particles and events. 
The first two harmonics vi and V2 arc referred to as di- 
rected and elliptic flow, respectively. 

It has been found that at low pT [pt ~ 2.0 GeV/c) 
the magnitude and trends of 2 are under-predicted by 
hadronic cascade models supplemented with string dy- 
namics [Tsj . but are well reproduced by models which 
either incorporate hydrodynamic flow [7|, |9| with a first 
order phase transition and rapid thermalization, r ~ 1 
fm/c [19|, or use a quasi-particle ansatz but include more 
than just 2-to-2 interactions [20l |. 

The mass dependence of V2 as a function of pt has been 
studied using identified baryons and mesons [l9l . [2l| and 
empirical scaling of elliptic flow per constituent quark 
was observed when the signal and the pT of the hadron 
were divided by the number of constituent quarks Uq 
(uq = 2 for mesons, 3 for baryons). This scaling is 
most clearly observed by plotting the data as a func- 
tion of transverse kinetic energy KEy = — m = 
\J Pt^ + 'Ti^ — m (23 | , where and to denote the trans- 
verse mass and mass of the particle, respectively. A re- 
cent study [2^ finds that the constituent quark scaling 
holds up to KEji w 1 GeV. This indicates partonic, 
rather than hadronic fiow, and suggests that the bulk 
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matter collectivity develops before hadronization takes 
place [H, m, [2^ . Results for the of the </> meson fur- 
ther validate the observation of partonic collectivity. The 
is not expected to be affected by hadronic interactions 
in the late stages of the medium evolution, due to its 
small interaction cross section with non-strange hadrons 

[la. 

All of the vi measurements referenced above were per- 
formed using the event plane method [1^. In PHENIX 
the event plane was determined at forward and back- 
ward pseudorapidities (jryl = 3.1-3.9) with the assump- 
tion that correlations induced by elliptic fiow dominate 
over all other non-fiow correlations [1^. Non-fiow corre- 
lations are those which are not correlated with the reac- 
tion plane. Common sources of non-flow correlations in- 
clude jets, the near-side ridge, quantum correlations and 
resonance decays. Simulation studies [l9l.[29| have shown 
that the correlations from jets and dijets become negli- 
gible when the rapidity separation between the particles 
and the event plane is greater than three units. Thus 
we expect that the event plane at forward pseudorapidi- 
ties I77I = 3.1-3.9 in the PHENIX experiment would not 
have significant jet-correlation with particles measured 
within the PHENIX central arm spectrometer covering 
the pseudorapidity window |?7| < 0.35. PHOBOS has ob- 
served that in central Au+Au collisions there is a ridge 
of particles (30j that are correlated in azimuthal angle 
with a high-pT particle and that this ridge of particles 
extends to I77I < 4 (for mid-rapidity triggers). The ridge 
could produce a non-fiow correlation that we can provide 
information by using our vi measurements that are made 
with different techniques and at different rapidities. 

Event-by-event flow fluctuations can also affect the 
magnitude of the extracted flow signal [3l|. This oc- 
curs because the event plane at forward pseudorapidities 
is reconstructed using particles from participant nuclcons 
whose positions fluctuate evcnt-by-event. Assuming that 
1)2 fluctuates according to a Gaussian distribution, the vi 
fluctuation is proportional to the fluctuation of the initial 
geometry. This effect scales as 1/A'^part, where A^part de- 
notes the number of participant nucleons. The difference 
between values obtained from different methods can 
be quantitatively understood in terms of non-flow and 
fluctuation effects [s^ . 

Hence in this paper we will compare the vi results from 
the event plane determined at two different pseudorapidi- 
ties with the goal to investigate the sensitivity of to 
non-flow and fluctuation effects. Additionally, we extract 
the elliptic flow with the two-particle cumulant method, 
which is expected to have higher sensitivity to non-flow 
contributions to Vi- 

In this paper, we describe the PHENIX measurements 
of elliptic flow (v2) at midrapidity (jTyj < 0.35) in Au-|-Au 
collisions at ^s„„ = 200 GeV obtained from a cumulant 
analysis of two-particle azimuthal correlations and the 
event plane method over a broad range of pt (pt = 0.2- 
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8 GeV/c) and ccntrality (0-60%). The paper is organized 
as follows: Section |TT] describes the PHENIX apparatus, 
with an emphasis on the detectors relevant to the pre- 
sented results, as well as the track selections used in the 
analysis. Section Hill gives details of the event-plane and 
the cumulant methods as applied in PHENIX. Section llVl 
discusses the systematic uncertainties of the event-plane 
and cumulant methods. The results from the two meth- 
ods are reported in Section |Vl Section |Vl] presents a 
comparison of V2 results across different experiments and 
discussion. The V2 values obtained from the different 
methods are tabulated in the Appendix. 



II. EXPERIMENTAL ANALYSIS 
A. The PHENIX detector 




West Beam View East 




FIG. 1: PHENIX experimental layout in 2004. The top panel 
shows the PHENIX central arm spectrometers viewed along 
the beam axis. The bottom panel shows a side view of the 
PHENIX muon arm spectrometers. 

The PHENIX detector consists of two central spec- 
trometer arms at midrapidity that are designated East 
and West for their location relative to the interaction re- 
gion, and two muon spectrometers at forward rapidity, 
similarly called North and South. A detailed descrip- 
tion of the PHENIX detector can be found in Ref. [33j . 



The layout of the PHENIX detector during data taking 
in 2004 is shown in Fig. [T] Each central spectrometer 
arm covers a pseudorapidity range of \rj\ < 0.35 subtend- 
ing 90 degrees in azimuth and is designed to detect elec- 
trons, photons and charged hadrons. Charged particles 
are tracked by drift chambers (DC) positioned between 
2.0 m and 2.4 m radially outward from the beam axis 
and layers of multi-wire proportional chambers with pad 
readout (two in the east arm and three in the west arm) 
PCI, PC2 and PC3 located at a radial distance of 2.4 
m, 4.2 m and 5 m, respectively. Particle identification 
is provided by Ring Imaging Cerenkov counters (RICH), 
a time-of- flight scintillator wall (TOP), and two types of 
electromagnetic calorimeters (EMCAL), the lead scintil- 
lator (PbSc) and lead glass (PbGl). 

The detectors used to characterize each event are the 
beam-beam counters (BBCs) [s^l and the zero-degree 
calorimeters (ZDCs) |35[. These detectors arc used to de- 
termine the time of the collision, the position of the colli- 
sion vertex along the beam axis and the collision central- 
ity and also provide the event trigger. In this analysis the 
BBCs are also used to determine the event plane. Each 
BBC is composed of 64 elements and a single BBC ele- 
ment consists of a one-inch diameter mesh dynode pho- 
tomultiplier tube (PMT) mounted on a 3 cm long quartz 
radiator. The BBCs are installed on the north and south 
sides of the collision point along the beam axis at a dis- 
tance of 144 cm from the center of the interaction region 
and surround the beam pipe. The BBC acceptance cov- 
ers the pseudorapidity range 3.1 < Ir]] < 3.9 and the full 
range of azimuthal angles. 

The ZDCs are hadronic calorimeters located on both 
sides of the PHENIX detector. Each ZDC is mechanically 
subdivided into 3 identical modules of two interaction 
lengths. They cover a pseudorapidity range of \r]\ > 6.5 
and measure the energy of the spectator neutrons with 
a 20 GeV energy resolution (35j . The shower maximum 
detectors (ZDC-SMDs) are scintillator strip hodoscopes 
between the first and second ZDC modules. This loca- 
tion approximately corresponds to the maximum of the 
hadronic shower. The horizontal coordinate is sampled 
by 7 scintillator strips of 15 mm width, while the vertical 
coordinate is sampled by 8 strips of 20 mm width. The 
active area of a ZDC-SMD is 105 mm x 110 mm (hori- 
zontal X vertical dimension). Scintillation light is deliv- 
ered to a multichannel PMT M16 by wavelength-shifter 
fibers. The ZDC-SMD position resolution depends on 
the energy deposited in the scintillator. It varies from < 
3 mm when the number of particles exceeds 100, to 10 
mm for a smaller number of particles. 



B. Event selection 

For the analyses presented here we used approxi- 
mately 850 X 10^ minimum-bias triggered events. The 
minimum-bias trigger was defined by a coincidence be- 
tween North and South BBC signals and an energy 
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FIG. 2: The correlation between ZDC energy and BBC 
charge sum for Au + An collisions at ysjy„ = 200 GeV. 
Solid lines represent the corresponding centrality boundaries 
up to 60% centrality bin. 



threshold of one neutron in the ZDCs. The events are 
selected offline to be within a z-vertex of less than 30 cm 
from the nominal center of the PHENIX spectrometer. 
This selection corresponds to 92.2~!il'^% of the 6.9 barn 
Au+Au inelastic cross section at = 200 GeV [36j . 

The event centrality was determined by correlating the 
charge detected in the BBCs with the energy measured 
in the ZDCs, as shown in Fig. [S] 

A Glauber model Monte-Carlo simulation [s^, [H, [s^ 
that includes the responses of BBC and ZDC gives an 
estimate of the average number of participating nucleons 
(iVpart) for each centrality class. The simulation did not 
include fluctuations in the positions of the nucleons which 
give rise to eccentricity fluctuations. Table |T] lists the 
calculated values of (A^part) for each centrality class. 



TABLE I: Centrality classes and average number of partici- 
pant nucleons (Afpart) obtained from a Glauber Monte-Carlo 
simulation of the BBC and ZDC responses for Au-|-Au colli- 
sion at y^sTnv = 200 GeV. Each centrality class is expressed 
as a percentage of (tauAu = 6.9 b inelastic cross section. Er- 
rors denote systematic uncertainties from the Glauber MC 
simulation. 



Centrality 



(A"part) 



0-10% 
10-20% 
20-30% 
30-40% 
40-50% 
50-60% 



325.2 ±3.3 
234.6 ±4.7 
166.6 ± 5.4 
114.2 ±4.4 

74.4 ±3.8 

45.5 ±3.3 



Track selection 



Charged particle tracks are measured using informa- 
tion from the DC, PCI and PC3 detectors and the z- 
vertex from the BBC. The DC has 12 wire planes which 
are spaced at 0.6 cm intervals along the radial direction 
from the beam axis. Each wire provides a track position 
measurement, with better than 150 fim spatial resolution 
in the azimuthal (0) direction. The PCI provides a space 
point in the (j) a-nd beam directions, albeit with lower res- 
olution. This space point and the vertex position help de- 
termine the three-dimensional momentum vector by pro- 
viding the polar angle for charged tracks at the exit of 
the DC. Trajectories are confirmed by requiring matching 
hits at PCS to reduce secondary background. Tracks are 
then projected back to the collision vertex th roug h the 
magnetic field to determine the momentum p [40| . The 
momentum resolution is Sp/p ~ 0.7%© 1.0% xp (GeV/c). 
The momentum scale is known to 0.7%, as determined 
from the reconstructed proton mass using the TOF de- 
tector. Further details on track reconstruction and mo- 
mentum determination can be found in Refs. [s^. [40l|. 

The tracks reconstructed by the DC which do not orig- 
inate from the event vertex have been investigated as 
potential background to the charged particle measure- 
ment. The main background sources include secondary 
particles from decays and e~^e~ pairs from the conver- 
sion of photons in the material between the vertex and 
the DC [1^. Tracks are required to have a hit in the 
PC3, as well as in the EMCAL, within at most 2 ct of the 
expected hit location in both the azimuthal and beam 
directions. This cut reduces the background not origi- 
nating in the direction of the vertex. In order to reduce 
the conversion background we further require tracks to 
have E/pt > 0.2, where E denotes the energy deposited 
in the EMCAL and px is the transverse momentum of 
particles measured in the DC. Since most of the electrons 
from photon conversion are genuine low pT particles that 
were reconstructed as high px particles, requiring a large 
deposit of energy in the EMCAL suppresses the electron 
background [4l|. We also require that there are no as- 
sociated hits in the RICH. The RICH is filled with CO2 
gas at atmospheric pressure and has a charged particle 
threshold 7th = 35 to emit Cerenkov photons. 



III. METHODS OF AZIMUTHAL ANISOTROPY 
MEASUREMENT 



In this section we introduce the methods for azimuthal 
anisotropy measurements as used in the PHENIX exper- 
iment. Section HI. A describes the event plane method 
using the BBCs and ZDC-SMDs detectors and Sec. III.B 
describes the two-particle cumulant method. 
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A. Event plane method 

The event plane method [1^ uses the azimuthal 
anisotropy signal to estimate the angle of the reaction 
plane. The estimated reaction plane is called the "event 
plane" and is determined for each harmonic of the Fourier 
expansion of the azimuthal distribution. The event flow 
vector Qn ~ iQx,Qy) and azimuth of the event plane 
for n-th harmonic of the azimuthal anisotropy can be 
expressed as 



M 



Qx = \Qn\cos{n'^n) ^^Wi cos (ncpi), 



M 



Qy = IQnl sin sin 
= — tan" 



1 fQy 



n 



(2) 

(3) 
(4) 



where M denotes the number of particles used to deter- 
mine the event plane, is the azimuthal angle of each 
particle, and Wi is the weight chosen to optimize the event 
plane resolution. Once the event plane is determined, the 
elliptic flow V2 can be extracted by correlating the az- 
imuthal angle of emitted particles with the event plane 



,,obs 



V2 



(cos (2 [0 -*„])) 



(5) 



Rcs{*„} (cos(2[^'„-*Rp]))' 

where (f> is the azimuthal angle of tracks in the laboratory 
frame, '^n is the n-th order event plane and the brackets 
denote an average over all charged tracks and events. The 
denominator Res{^'„} is the event plane resolution that 
corrects for the difference between the estimated event 
plane and true reaction plane ^'rp. 

In this paper the second-harmonic event planes were 
independently determined with two BBCs located at for- 
ward (BBC South, referred to as BBCS) and backward 
(BBC North, referred to as BBCN) pseudorapidities \rj\ 
= 3.1-3.9 [19[. The difference between the two indepen- 
dent event planes was used to estimate the event plane 
resolution. The planes were also combined to determine 
the event plane for the full event. A large pseudorapidity 
gap between the charged particles detected in the cen- 
tral arms and the event plane at the BBCs reduces the 
effect of possible non-flow contributions, especially those 
from dijets [2^. The measured V2 of hadrons in the cen- 
tral arms with respect to the combined second-harmonic 
BBC event plane will be denoted throughout this paper 
as wzlBBC}. 

Two first-harmonic event planes were also determined 
using spectator neutrons at the two shower maximum 
detectors (ZDC-SMDs) that are sandwiched between the 
first and second modules of the ZDCs. Forward (ZDCS) 
and backward (ZDCN) SMDs which cover pseudorapid- 
ity |?7| > 6.5 were used. The measured V2 of hadrons in 
the central arms determined with respect to the first- 
harmonic ZDC-SMD event plane will be denoted as 
i;2{ZDC-SMD}. 



The pseudorapidity gap between the hadrons measured 
in the central arms and the ZDC-SMDs is larger than 
that for the BBCs which could cause a further reduction 
of non-flow contributions on z;2{ZDC-SMD}. Since the 
ZDC-SMD measures spectator neutrons, the ZDC-SMD 
event plane should be insensitive to fluctuations in the 
participant event plane. Hence fluctuations in i'2{ZDC- 
SMD} should be suppressed up to fluctuations in the 
spectator positions. 

For completeness, two further event planes are defined 
1) a combined event plane defined by the weighted aver- 
age of event planes at the forward and backward pseu- 
dorapidities for both BBCs and ZDC-SMDs, and 2) an 
event plane found using tracks in the central arm. The 
event plane at the central arms (CNT) is only used to es- 
timate the resolution of BBC and ZDC-SMD event planes 
by using three subevents combination of the ZDC-SMD, 
BBC and CNT. 



1. Event plane determination 

To determine an event plane the contribution at each 
azimuthal angle needs to be appropriately weighted de- 
pending on the detector used. For the BBC we chose 
the weights to be the number of particles detected in 
each phototube, while for the ZDC-SMD the weights were 
based on the energy deposited in each of the SMD strips. 
For the CNT event plane the weight was taken to be 
proportional to pt up to 2 GeV/c and constant for pT > 
2 GeV/c. For the CNT event plane we also adopted a unit 
weight {wi = 1) and found that the resulting CNT event 
plane resolution extracted by comparing the CNT event 
plane with the BBC and ZDC-SMD planes was nearly 
identical when using the pr-dependent or unit weights. 

Corrections were performed to remove possible biases 
from the finite acceptance of the BBC and ZDC-SMD. 
In this analysis, we applied two corrections called the re- 
centering and shift methods. In the re-centering method, 
event flow vectors are shifted and normalized to a Gaus- 
sian distribution by using the mean (Q) and width a of 
flow vectors; 



Qx {Qx} 



Qy {Qy) 



(6) 



This correction reduces the dependence of the event plane 
resolution on the laboratory angle. Most acceptance ef- 
fects were removed by the application of the re-centering 
method. However, remaining small corrections were ap- 
plied after re-centering using the shift method [28| . in 
which the reaction plane is shifted by A^Pn defined by 

nA*„(«'„) = ^ -[- (sin(fcn*„))cos(fcn*„) 

k=l 

+ (cos (fcn*„)) sin (fcn*„)], (7) 

where fc^ax = 8 in this analysis. The shift ensures 
that dN/d^n is isotropic. When k^sx was reduced to 
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fcmax = 4, the difference in the extracted V2 was negligi- 
ble and thus we include no systematic uncertainty due to 
the choice of fcmax in our V2 results. 

Independent corrections were applied to each centrality 
selection in 5% increments and in 20 cm steps in z-vertex 
in order to optimize the event plane resolution. The cor- 
rections were also done for each experimental run (the 
duration of a run is typically 1-3 hours) to minimize the 
possible time-dependent response of detectors. 
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FIG. 3: Event plane distributions after applying all correc- 
tions for the ZDC-SMD (circles), BBC (triangles) and CNT 
(squares). The statistical error bars are smaller than the sym- 
bols. The distributions for the BBC and CNT event planes 
are scaled by 3/4 and 1/2 to improve visibility. 

Figure [3] shows event plane distributions for a sub- 
sample of the entire data set. After all corrections are 
applied the event plane distributions are isotropic. 



2. Event plane resolution 

The event plane resolution for V2 was evaluated by both 
the two-subevents and three-subevents methods. In the 
two-subevents method the event plane resolution [2^ is 
expressed as 



(cos(A;n[«'„ - *rp])) = 



2^/2 



Xne 



I(k+l)/2 



(8) 



where Xn = VnV^M, M is the number of particles used 
to determine the event plane 5'„, Ik is the modified Bessel 
function of the first kind and k = 1 for the second har- 
monic BBC event plane. For the ZDC-SMD event plane 
the resolution is estimated with both fc = 1 or 2 in Eq. ([S]). 
We will discuss the difference between these estimates in 
Sec HvH 

To determine the event plane resolution we need to 
determine Xn- Since the North and South BBCs have 



approximately the same rj coverage, the event plane res- 
olution of each sub-detector is expected to be the same. 
The same is true for the North and South ZDC-SMDs. 
Thus, the subevent resolution for South and North event 
planes can be expressed as 



cos 



(2[VI/S(N) „ vl^j^p])^ = y/(eOS (2[VI/S _ VI/N])), (9) 



where ^'n*'^'' denotes the event plane determined by the 
South (North) BBC or ZDC-SMD. Once the subevent 
resolution is obtained from Eq. one can calculate 
X^n^ using Eq. The Xn for the full event can then 
be estimated by Xn = V^Xn"*^- This is then substituted 
into Eq. ([5]) to give the full event resolution. Since the 
multiplicity of the full event is twice as large as that of 
the subevent, Xn is proportional to ^/M. 

In the three-subevents method the resolution of each 
subevent is calculated by adding a reference event plane 
^-^^ in Eq. dl]): 



Res{vl/f} = J(cos(2[vI/A_vl.Bj)) 



^ (cos(2[v&g-^A])^ 
(cos(2[vI/B _VI/C]))' 



(10) 



where 1, 111,71 are the harmonics of the event plane for 
subevent A, B and C, respectively. The multiplicity of 
each subevent is not necessarily the same in Eq. (fTUl) . 

The resolution of each sub-detector for the BBC and 
ZDC-SMD can be evaluated with the three-subevents 
method. For the BBC event plane the reference event 
plane is chosen to be the ZDC-SMD event plane and vice 
versa. We found that the agreement of the event plane 
resolutions for BBCS and BBCN is much better than 1%, 
while the ZDCS and ZDCN resolutions are comparable 
with each other within 2%. 

Figure 0] shows the full-event resolution as a func- 
tion of centrality. The resolution of ZDC-SMD is much 
smaller than that of BBC because the resolution of the 
first-harmonic event plane is proportional to (xi)^- The 
dashed lines are the resolutions obtained from the three- 
subevents method with the CNT event plane as the refer- 
ence plane. For example, the BBC event plane resolution 
is estimated by substituting *f *f 



and *C ^ ^ZDC-SMD gy including the 

CNT event plane, the BBC resolution increases by about 
3% compared to that of the two-subevents method. For 
the ZDC-SMD we observe the opposite effect, namely 
the resolution decreases by about 8%. In Sec. IVII the re- 
suhing U2{BBC} and U2{ZDC-SMD}, corrected by the 
resolution obtained using the ZDC-BBC-CNT combina- 
tion, will be compared to those with the resolution deter- 
mined from South-North subevents. Table HIl summarizes 
the event plane resolutions. 
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FIG. 4: Full-event resolutions for the ZDC-SMD (filled cir- 
cles) and BBC (open diamonds) from the two-subevents 
method, Eq. ((8)1, as a function of centrality in Au-|-Au at 
^s„„ = 200 GeV. The dashed lines represent resolutions 
from the three-subevents method with the CNT event plane 
as a reference. Statistical errors are smaller than the symbols. 



TABLE II: Event plane resolutions for centrality 0-60% at 
^/sj^ = 200 GeV. S-N denotes the resolutions estimated 
from South and North correlation of BBC and ZDC-SMD us- 
ing Eq. (O and dU, and resolutions for ZDC-BBC-CNT are 
estimated from Eq. (|10|) . The errors are statistical only. 



Res{*?BC| 



Centrality 



S-N 



ZDC-BBC-CNT 



0-10% 
10-20% 
20-30% 
30-40% 
40-50% 
50-60% 


0.2637 ± 0.0003 
0.3809 ± 0.0002 
0.3990 ± 0.0002 
0.3634 ± 0.0002 
0.2943 ± 0.0003 
0.2106 ± 0.0004 


0.272 ± 0.003 
0.394 ± 0.001 
0.4106 ± 0.0008 
0.3759 ± 0.0007 
0.3067 ± 0.0007 
0.2240 ± 0.0009 


Res{v[>™c-SMDj 


Centrality 


S-N 


ZDC-BBC-CNT 


0-10% 
10-20% 
20-30% 
30-40% 
40-50% 
50-60% 


0.02 ± 0.01 
0.059 ± 0.003 
0.087 ± 0.002 
0.100 ± 0.002 
0.102 ± 0.002 
0.100 ± 0.002 


0.0223 ± 0.0003 
0.0574 ± 0.0002 
0.0818 ± 0.0002 
0.0928 ± 0.0002 
0.0920 ± 0.0002 
0.0798 ± 0.0003 



3. Correlation of event planes 

Figure O shows the correlation of two different event 
planes as a function of centrality. The first harmonic 
event plane correlation for South-North detector com- 
binations is negative both for the ZDC-SMDs and the 
BBCs over all centrality bins, as shown in Fig. [^a). 
This is due to the fact that vi is an odd function of rj. 
The magnitude of the ZDC-SMDs correlation is about a 
factor of two larger than that of the BBCs for midccn- 
tral collisions. This indicates that the magnitude of vi 
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FIG. 5: (a) Correlation of first harmonic event planes between 
forward and backward ZDC-SMDs (filled circles) and BBCs 
(open diamonds) as a function of centrality. (b) Correlation 
of first harmonic event planes between ZDC-SMDs and BBCs 
as a function of centrality, where filled (open) squares are 
the correlation for opposite side (same side) of 77 subevents. 
Statistical errors are smaller than the symbols. 



and/or the subevent multiplicity at higher pseudorapidi- 
ties are larger compared to that at the BBC location, 
since the magnitude of the correlation is proportional to 
ViAI. Fig. Ob) shows the correlation of the first har- 
monic event planes between BBC and ZDC-SMD. The 
same-side rj correlation is negative while the opposite-side 
r] correlation is positive, which shows that the particles 
detected at the BBCs (dominantly charged pious emit- 
ted from participant nucleons) have the opposite sign of 
vi compared to the spectator neutrons detected at the 
ZDCs-SMDs. 

The correlation of the mixed harmonic event planes 
provides the sign of V2 since the correlation is given by 
the expression [23 



(C0S(2[VI/ZDC-SMD_^BBC])^ 



- (Res{vI/fDC-SMD|)^ j^gg|^BBC| 



±2V2- (cos(*? 



ZDCS 



ZDCN 



)) 



X y^(c0s(2[*fBCS_^BBCN])^_ (i^) 



Three assumptions were made to obtain Eq. (fTT|) : (1) the 
BBC and ZDC-SMD are statistically independent, (2) 
the weak flow limit is applicable, and (3) the subevent 
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FIG. 6: The correlation between the first harmonic ZDC-SMD 
and the second harmonic BBC event planes as a function of 
centrality. The dashed line shows the result obtained using 
Eq. (|lip . Statistical errors are smaller than the data symbols. 



multiplicity Al is equal in the North-South direction for 
the same detector type. Thus the sign of the correlation 
of the mixed harmonic event planes in Eq. (|lip is deter- 
mined by the term Res{\I'2 ^^}, which in turn determines 
the sign of V2 measured at the BBC. 

Figure [6] shows the mixed harmonic correlation of the 
ZDC-SMD and BBC event planes as a function of cen- 
trality. The approximations in Eq. (|lip provide a good 
description of the magnitude of the measured correlation 
as shown by the dashed line. The correlation is positive 
over all centrality bins. This result indicates that the 
sign of V2 at the BBC is positive. 



B. Cumulant method 

In this section, we present the application of the cu- 
mulant method for azimuthal anisotropy measurements 
in PHENIX. This method uses cumulants of multiparticle 
correlations [l^, to extract the azimuthal anisotropy. 
The cumulant method has been successfully applied in 
several heavy-ion experiments utilizing detectors with 
full azimuthal coverage (NA49, STAR) HSll. Here, we 
describe the first application of the method for a detec- 
tor with only partial azimuthal coverage. The cumulant 
method does not require the measurement of the reaction 
plane, instead the cumulants of multi-particle azimuthal 
correlations are related to the flow harmonics u„, where 
n is the harmonic being evaluated. The cumulants can 
be constructed in increasing order according to the num- 
ber of particles that are correlated with each other. Since 
PHENIX has partial azimuthal coverage, reliable extrac- 
tion of azimuthal anisotropy requires the choice of a fixed 
number of particles from each event in order to avoid ad- 
ditional numerical errors [i^ . 

Particles in an event are selected over a fixed pT range 



where there is sufficient multiplicity. These particles 
(called "integral particles" hereafter) are used to deter- 
mine integrated flow, that is flow measured over a large 
{pt,ti) bin. For differential flow measurement, we select 
particles (called "differential" particles) over small {pt,ii) 
bins, from which the integral particles are excluded so as 
to avoid autocorrelations. For each event a fixed number 
M of particles, chosen at random among the integral par- 
ticles in the event, are used to reconstruct the integrated 
flow through the generating function G2{z) defined by: 

M 

where Wj is the weight, chosen to be equal to 1 in our 
analysis, 4>j is the azimuth of the detected particles, and 
M is the multiplicity chosen for the integrated flow re- 
construction. 6*2(2:) is a function of the complex variable 
z. The average of 6*2(2) over events is then expanded in a 
power series to generate multi-particle azimuthal correla- 
tions. The generating function of the cumulants, defined 
by 



C„(z) = M 



((G„(z))' 



/M 



- 1 



(13) 



generates cumulants of azimuthal correlations to all or- 
ders, the lowest being the second order, as detailed in 
Section II. B of Ref. [4^. The formulas used to compute 
the cumulants from which the V2 is computed are given 
in Appendix B of Ref. [i^. In the case of a perfect ac- 
ceptance the relations between the anisotropy parameter 
V2 and the lowest order cumulants are 



V2{2Y = C2{2}, 
V2{AY = -C2{4}, 



(14) 
(15) 



for the integrated anisotropy. Here W2{2} and W2{4} are 
the second and fourth order 112, respectively; whereas, 
C2{2} and C2{4} are the second and fourth order cumu- 
lants. Because the typical multiplicity of charged hadrons 
in PHENIX did not allow a reliable calculation of W2{4}, 
we report here only the 1^2 {2} results. 

The remaining differential particles in the same event 
are selected in different {pt, v) bins and the differential 
cumulants are calculated from the generating function 



V2/2iz) 



^ {e^''^G2{z)) 
~ {G2{z)) ■ 



(16) 



where {G{z)) denotes an average over all events, and ip is 
the azimuth of each differential particle. D2/2 denotes the 
second order differential cumulant computed with respect 
to the second order integral cumulant. 

The differential W2/2{2}(pt, ??), the second order differ- 
ential 2 with respect to the second order integrated 2, is 
calculated from the relation 



V2/2{'2}(PT,V) 



d2/2{2}{pT,v) 
V2{2} ■ 



(17) 
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where d2/2{'2}{PT,v) is the second order differential cu- 
mulant. These relations have to be modified through 
acceptance corrections which are detailed below. 



1. Acceptance/efficiency corrections 

The central arms detectors in PHENIX have only par- 
tial azimuthal coverage and the implementation of the 
cumulant method requires an additional acceptance cor- 
rection. In order to correct for the influence of the de- 
tector acceptance on the raw anisotropy values, we ap- 
ply a correction factor using the prescription described in 
Ref. [13 . The acceptance and efficiency of the detector is 
characterized by a function pT, rj) which is expressed 
in terms of the Fourier series 



+ 00 

p— — oo 



(18) 



The Fourier cocfhcients ap{pT,'ri) for the detector accep- 
tance were extracted from the fit of the respective az- 
imuthal distributions of integral and differential parti- 
cles. The coefficients resulting from such fits were then 
used to calculate the correction factor for the raw values 
of the V2 following the procedure detailed in Appendix C 
of Ref. m. 



centrality = 20 - 40% 
1.2<p.|.<1.4 GeV/c 

Fourier fit 




300 

(degrees) 



FIG. 7: Azimuthal angular distribution and corresponding 
Fourier fit for centrality 20-40% and pr = 1.2-1.4 GeV/c. 

Figure[7]shows a typical azimuthal angular distribution 
of differential particles detected in the PHENIX central 
arms and the corresponding Fourier fit used to correct for 
acceptance inhomogeneities. The Fourier fit reproduces 
well the overall features of the acceptance profile. This 
produces typical correction factors that are in the range 
1.1-1.2 for the differential flow and depend very little on 
centrality and pr, as shown in Fig. [S] 
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FIG. 8: (a) Acceptance correction factor for differential f2{2} 
as a function of pr for centrality 10-20% (b) Acceptance cor- 
rection factor as a function of centrality for pr range 0.4-0.5 
GeV/c in Au+Au collisions at ^/s^ = 200 GeV. 



2. Simulations 

While Fig. [7] shows that the uneven detector accep- 
tance is reproduced by the Fourier fit, a better test of 
the cumulant method is to use Monte-Carlo simulations, 
as was done in Ref. [1^ . For these tests events were gen- 
erated with particles having a distribution of the form 
1 -I- 2vi cos (j) + 2v2 cos 20, with known integrated and dif- 
ferential azimuthal anisotropics. The anisotropy was in- 
troduced into the events by way of a Fourier weighted 
selection of the azimuthal angles followed by a random 
event rotation designed to simulate the random orien- 
tation of the reaction plane. The multiplicity of these 
events was chosen to reflect the typical multiplicity mea- 
sured with the PHENIX detector and the angles were 
chosen from a fllter that is representative of the PHENIX 
acceptance. We extracted Fourier components from these 
simulated results and applied these to extract corrected 
elliptic flow values. 

Figure [H] shows selected results from these simulations. 
Corrected differential anisotropy values are compared for 
various input differential V2 values, with the integral V2 
kept fixed. The dotted line shows the trend expected if 
the extracted V2 is identical to the input value used to 
generate the events. The good agreement between the 
input and extracted V2 attests to the reliability of the 
analysis method within the acceptance of the PHENIX 
central arms. 



IV. SYSTEMATIC UNCERTAINTIES 

In this section, we present the systematic uncertainties 
on the V2 from the event plane method (Sec. IIV Ap and 
the two-particle cumulant method (Sec. IIVB|) . Table Hill 
lists the different sources of systematic errors for each 



11 




0.04 0.06 
V2(input) 



0.1 



FIG. 9: Comparison of input and extracted differential V2 val- 
ues for a fixed integral V2 of 8%. The dotted line indicates the 
expectation if input and reconstructed values are the same. 
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FIG. 10; The radial PCS matching distribution for real (open 
circles) and random tracks (solid lines) for 6 < pr < 8 GeV/c 
in centrality 0-60%. 



method. The errors in Tab. IIIII are categorized by type: 

A point-to-point error uncorrelated between px bins, 

B pt correlated. aU points move in the same direction 
but not by the same factor, 

C an overaU normalization error in which all points 
move by the same factor independent of px- 



TABLE III: List of systematic uncertainties given in percent 
on the ii2{ZDC-SMD}, W2{BBC} and f2{2} measurements. 
The ranges correspond to different systematic errors for dif- 
ferent centrality bins. 



Error source 


Percentage error 


Type 


t;2{BBC} i;2{ZDC-SMD} 


Background contribution 


< 5% in PT < 4 GeV/c 


B 




5-30% in PT > 4 GeV/c 


B 


Event plane calibration 


1-5% 


C 


Event plane determination 


1-4% 1-16% 


C 


Acceptance effect 


1% 1-25% 


c 


on event planes 






V2{2} 


Fixed multiplicity 


5% 


B 


Integrated pr range 


3-8% 


B 


Background correction 


6-10% 


B 



A. Event plane method 

1. Background contributions 

In order to study the influence of background on our 
results, we varied one of the track selections while keeping 
other cuts fixed and investigated the effect on V2 in the 
following two cases: (i) the PC3 and EMCAL matching 
cuts, ± 1.5 and ± 2.5a matching cuts and (ii) E > QAbpr 



and E > 0.25pT- For both conditions, we found that the 
difference of the V2 is 1-2% for pt < 4 GeV/c, and 5-20% 
for Pt > 4: GeV/c depending on px and centrality. 

The effect of the RICH veto cut has also been studied. 
Since the contribution of charged tt increases without the 
RICH veto cut, the p/n ratio decreases at highpT- Thus, 
the V2 for charged hadrons could be modified due to the 
difference of f 2 between protons and tt in the range 4 < 
PT < 8 GeV/c. We found that V2 is 10-20% different 
without the RICH veto cut for pr > 4-5 GeV/c, where 
the charged tt starts firing the RICH. 

One of the remaining sources of background contribu- 
tion comes from the random tracks that are accidentally 
associated with the tracks in PC3. These random tracks 
have been estimated by swapping the z-coordinate of the 
PC3 hits and then by associating those hits with the real 
tracks. Figure [TO] shows the comparison of the radial 
PCS matching distribution between the real and random 
tracks for 6 < pt < 8 GeV/c. The signal to background 
ratio S/ B is evaluated in the crpc3 < 2 window, and is ^ 
52 for 6 < PT < 8 GeV/c in centrahty 0-60%. 

The ratio of real and random tracks with and without 
the E/pT > 0.2 cut is shown as a function of pr for 
centrality 0-60% in Fig. [Ill The E/pr > 0.2 cut reduces 
the random tracks and improves the S/ B ratio by a factor 
of w 10-24 for Pt > 4 GeV/c. Since random tracks are 
not expected to be correlated with the event plane, we 
assume that their V2 ~ and evaluate the systematic 
uncertainty on V2 to be less than 2% for pT > 4 GeV/c, 
increasing to 5% for pT < 0.5 GeV/c. 

There is a finite residual background contribution even 
after the E/pT > 0.2 has been applied, as observed in 
Fig. [TOl The residual backgrounds have been estimated 
by fitting the apcs with a double Gaussian while requir- 
ing that the signal and residual background crpcs dis- 
tribution have the same mean. For the highest pt bin, 
we found that the signal to background ratio is ~ 5 for 
fpc3 < 2. The systematic error on V2 is evaluated by 
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FIG. 11: The ratio of real S to random tracks B as a function 
of pt in centrality 0-60%. Solid and open circles show the 
S/B ratio with and without E/pr > 0.2, respectively. 



comparing the measured V2 with that of signal 

B\ B 
'S 



(19) 



where t^f , V2 and V2 are respectively V2 of signal, back- 
ground estimated for crpcs > 3, and measured within 
the 2(7 matching window. The systematic uncertainties 
are less than 5% for < 4 GeV/c, and ^ 5-10% for 
higher pT- All the above systematic errors are added 
in quadrature and the overall systematic error from the 
background contribution is estimated to vary from < 5% 
for Pt < 4 GcV/c to --30% for higher pT- 



2. Event plane calibrations 

The procedures used in the determination and cali- 
bration of event planes are the dominant sources of sys- 
tematic errors on V2 and are discussed in the following 
sections. 

Different calibration procedures of the BBC event 
plane were extensively studied for previous Au -I- Au 
data sets [l^. We followed the same procedure to study 
the systematic errors on the BBC and ZDC-SMD event 
planes. Systematic uncertainties from the shift methods 
on 'i;2{BBC} are ^ 1-5% depending on the centrality. 
The systematic errors on the i;2{ZDC-SMD} are 1-2% 
larger than those on W2{BBC} for centrality 10-30% and 
50-60%, although those are still less than 5%. 



3. Event plane determination 

Figure [T^] shows the comparison of {V2) for different 
sub-detectors with respect to the BBC and ZDC-SMD 
event planes as a function of centrality. Systematic errors 



are estimated by taking the maximum difference of the 
V2 from the South and North event planes to that from 
the combined South- North event plane scaled by 2/ \/V2 
for each centrality. Systematic errors range from 1-4% 
for the BBC, and 1-16% for the ZDC-SMD event planes 
depending on the centrality bins. 



^. Effect of non-uniform acceptance on V2 



In this subsection we discuss the effect of non-uniform 
acceptance on the measured V2. In practice, the imper- 
fect azimuthal acceptance of the BBC or ZDC-SMD or 
the central arms could induce an azimuthal-dependent 
event plane resolution and/or smear the magnitude of 
V2. In order to study the possible effect of non- uniform 
acceptance, the measured V2 is decomposed into X and 
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FIG. 12: (a) Comparison of (U2) averaged over 0.2 < pt < 8 
GeV/c as a function of centrality for the BBC event planes. 
Open triangles and crosses represent the V2 with respect to the 
event planes from South and North sub-detectors and filled 
circles show the V2 from combined South-North event planes. 
Results from South and North event planes are shifted in the 
x-direction to improve visibility, (b) The same comparison 
for the ZDC-SMD event planes. Only statistical errors are 
shown and they are smaller than the symbols. 
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Y components [46 



V2 (cos(20)cos(2vl/^)) 
a+ Res{vl/A;X} ' 
V2(sin(2(/)) sin (2^;^)) 
Res{vl/A;Y} ' 



(20) 



where </> denotes the azimuthal angle of hadrons mea- 
sured in the central arms and af = l± (cos (4(^)) are 
the acceptance correction factors of the measured V2 in 
the central arms. The coefficient af should be unity in 
the case of perfect azimuthal acceptance. Rcs{5'^;X} 
and Res{\E'^;Y} denote the event plane resolution for 
and respectively and are expressed as 



Res{vI'f;X} = y'(cos(2vI/A)cos(2^FB))) 

/(cos (2^-0) cos (2*f)) 

V (cos(2vI/B)cos(2vI/C))' 

Res{vI'f;Y} = y^(sin(2vI/A)sin(2vI'B))) 

/ (sin(2vI/g)sin(2M/f)) 

V (Sin(2vl/B)sin(2vl/C))' 



(21) 



where l,m,n are the harmonics of event planes for 
subevents A, B, and C, respectively. Another accep- 
tance effect from the difference between Res{^'^; X} and 
Rcs{5'^; Y} is discussed below. 
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FIG. 13: Acceptance correction factors in the central arms 
as a function of pr for centrality 0-60%. Correction factors 
become unity for a perfect azimuthal acceptance. Statistical 
errors are smaller than the symbols. 



Figure fT3l shows the acceptance correction factor af as 
a function of pr in the central arms for centrality 0-60%. 
The pt dependence is parameterized by 



P2 



''(Pt) = 1 T Poe 



1 ^ e(PT~P3)/P<i 



P5 ,(22) 



where p„ (n = 0,1,. ..,5) are free parameters. From the fit, 
we get po = 0.131, pi = 1.203, p2 = 0.029, pz = 0.640, 
PA — 0.096 and p5 = —0.097. There is no centrality 
dependence of the acceptance corrections in the measured 
centrality range and these same correction factors are 
applied for all centrality bins. 
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FIG. 14: (a) Raw tJ2{BBC} without the acceptance correction 
as a function of pr in centrality 20-60% for vf (filled circles), 
(filled squares) with the South BBC event plane and or 
V2 (open circles), V2 (open squares) with the North BBC 
event plane, (b) The same comparison with the acceptance 
correction. 



Figure [14] shows the raw U2{BBC} as a function of 
Pt in 20-60% centrality bin. is systematically higher 
than for px > 1 GeV/c as shown in Fig. [TW a). Fig- 
ure [TH^b) shows that and agree with each other af- 
ter dividing Wj^* by , the remaining difference between 
them being accounted for as a systematic error. For the 
ZDC-SMD event plane we observed a similar trend for 
and uj. 

A possible non-uniform acceptance of the BBC 
and ZDC-SMD could lead to the difference between 
Res{5'„;X} and Res{^'„;Y}. If the azimuthal coverage 
of both detectors is perfect, Res{^'„; X} and Res{\I'„; Y} 
should be identical. Therefore, the effect of the accep- 
tance of the detector on the event plane resolution can 
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be assessed by comparing Res{^'„;X} and Res{^'„; Y}. 
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FIG. 15: (a) Comparison of Res{^'„;X} (open squares) and 
Res{*I'n;Y} (filled circles) with Res{>I',i} (dashed lines) for 
the BBC event plane (n = 2) as a function of centrality. The 
resolutions are calculated by using Eq. (|2ip with the ZDC- 
SMD, BBC and CNT event planes. Resj^n} is divided by 
\/2 in order to compare Res{»I'„;X} and Res{vl'„;Y}. (b) 
The same comparison for the ZDC-SMD event plane (n — 
1). Only statistical errors are shown and are smaller than 
symbols. 

Figure fTSl shows Res{*„;X} and Res{'I'„;Y} of the 
BBC and ZDC-SMD as a function of ccntrahty. The 
resolutions are calculated by using Eq. (PT|) with the 
ZDC-SMD, BBC and CNT event planes. Rcs{*„;X} 
was comparable with Res{^'„; Y} for both the BBC and 
ZDC-SMD event planes. They also agreed, within sta- 
tistical errors, with the expected resolution, namely the 
full event resolution scaled by 1/ a/2. We also evaluated 
Resj^-niX} and Resj^-n; Y} of BBC and ZDC-SMD for 
the two-subcvents method. Resj^'f X} was consis- 
tent with Rcs{*fi3C. Y}. However, for the ZDC-SMD 
event plane, Res{vI/ZDC-SMD. y} (Rcs{vI/ZDC-smd. ^j) 
was systematically higher (lower) by about 30% than 
the expected resolution when the resolutions were cal- 
culated with fc 1 in Eq. ([5]). The difference be- 
tween Res{vI/™c-SMD. ^nd Res{vI/ZDC-SMD. y} for 

the two-subevents method is attributed to the non- 
uniform acceptance between horizontal (x) and ver- 
tical (y) directions of the ZDC-SMD. Those reso- 
lutions of the ZDC-SMD were consistent with each 
other using k = 2. For k = 2, the non-uniform 
acceptance in the azimuthal directions cancels out 



since Res{«'i°^-^^°;X,Y} contain both (cos (*)) and 
(sin(^')) terms. Thus, Res{*i°^-™°; X, Y} should be 
the same and consistent with that from the expected res- 
olution. 
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FIG. 16: (a) Comparison of (open diamonds) and vj 
(open crosses) with the total 112 (filled circles) for the BBC 
event plane as a function of pr for the centrality bin 20-60%. 
Res{*I'„;X} and Res{'I'„;Y} are calculated by the combina- 
tion of the ZDC-SMD, BBC and CNT event planes. Accep- 
tance corrections are included into and . Error bars 
denote statistical errors, (b) The same comparison for the 
ZDC-SMD event plane. 



The comparison of and with V2 with respect 
to the BBC and ZDC-SMD event planes is shown in 
Fig. 1161 The maximum difference of and vj rela- 
tive to W2{BBC} is about 2% for the centrality range 
20-60% and is independent of centrality. Systematic un- 
certainties are evaluated by scaling the maximum differ- 
ence by 2/-\/l2. The same comparison is also made for 
W2{ZDC-SMD} as shown in the bottom panel in Fig. [161 
The systematic errors range from 1-25% and in this case 
strongly depend on the centrality, as well as on the cor- 
rections by the different event plane resolutions, and 
vj are ~ 10-25% different from i;2{ZDC-SMD} in the 0- 
20% centrality bin due to the very low resolution. This 
systematic uncertainty is denoted as "Acceptance effect 
on event planes" in Table WUl 
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B. Cumulant method 

The potential sources of systematic errors on the cu- 
mulant measurements are detailed below. 



1. Fixed multiplicity cut 

Following Ref. [i^ a fixed multiplicity is used to recon- 
struct the integrated flow to avoid introducing additional 
errors arising from a fluctuating multiplicity. In our anal- 
ysis the systematic errors were estimated by varying the 
fixed multiplicity cut used for the reconstruction of the 
integrated flow and studying its effect on the differential 
flow values. 
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FIG. 17: (a) V2{2} as a function of pt for centrality 10-20% 
in Au-|-Au collisions at ys„„ = 200 GeV for different fixed 
multiplicity cuts, corresponding to 60% (filled triangles) , 70% 
(open circles) and 80% (open crosses) of the mean multiplicity, 
(b) The ratio of V2{pt) for the two lowest multiplicity cuts to 
V2{pt) for 80% of the mean multiplicity. 



Figure [PTT a) shows the variation of V2 with pT for in- 
tegral multiplicity cuts equal to 60%, 70%, and 80% of 
the mean multiplicity for the centrality bin 20-40%. The 
ratio of the differential V2 values, shown in Fig. [TTT b). 
is used to estimate the systematic error on our measure- 
ments, which is ~ 5%. 



2. Pt range for integrated flow 

In order to assess the influence of the px range used to 
estimate the integrated flow on the differential flow, we 
chose different pT ranges over which the integral particles 
were selected. Differential V2 results were obtained for 
three px ranges: 0.25 - 2.0 GeV/c, 0.25 - 1.5 GeV/c and 
0.3 - 1.5 GeV/c. The systematic error from this source 
is estimated to be 3-8% depending on centrality and pt- 

3. Background contribution 

The procedures followed for studying the background 
contribution to f 2{2} were the same as for the event plane 
method. After background subtraction the systematic er- 
ror is calculated by determining the difference between 
the V2 obtained from using 2a and 3(t association cuts. 
We determined that the overall systematic error due to 
these differences is 6-10% depending on px and central- 
ity. 

V. RESULTS 
A. Pt dependence of V2 

The Pt dependence of V2 has been instrumental in 
revealing the hydrodynamic properties of the matter 
formed at RHIC 0, HTj. In this context, it is impor- 
tant to compare the px dependence of V2 from different 
methods to establish the robustness of our V2 measure- 
ments. This comparison is displayed in Fig. [TS] which 
shows the differential charged hadron V2 as a function of 
Pt from the event plane and cumulant methods for dif- 
ferent centrality bins in the range 0-60% in Au-|-Au at 
y/s^ = 200 GeV. f2{2} increases up to pT ~ 3 GeV/c 
and saturates at 0.1-0.25, depending on centrality, for 
higher pT- On the other hand, i'2{BBC} and W2{ZDC- 
SMD} reach their maximum value at pT ~ S GcV/c, and 
decrease for higher pT- 

The differences between -!;2{BBC} and W2{ZDC-SMD} 
are independent of pt within systematic errors in the 
measured centrality range. t>2{ZDC-SMD} is consistent 
with 'i;2{BBC} within systematic errors in the 0-40% cen- 
trality range, but is ^ 10-20% smaller than W2{BBC} in 
the 40-60% centrality range. These results could indi- 
cate that the influence of non-flow effects on D2{BBC} 
is small and within the systematic errors, because non- 
flow effects are not expected to influence W2{ZDC-SMD}. 
The difference between ?;2{BBC} and i;2{ZDC-SMD} in 
peripheral collisions could be attributed to non-flow con- 
tributions that might be proportionally larger in more 
peripheral collisions. 

The cumulant and event plane V2 agree well within 
systematic uncertainties in the centrality range 0-40%. 
In more peripheral collisions, there may be some differ- 
ences developing above pT — 4 GeV/c. Correlations be- 
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FIG. 18: Charged hadron V2{pt) in Au+Au collisions at ^/sj^ ~ 200 GeV from the two-particle cumulant method (filled 
squares), the BBC event plane (filled triangles) and the ZDC-SMD event plane (filled circles) for centrality (a) 0-10%, (b) 
10-20%, (c) 20-30%, (d) 30-40%, (e) 40-50%, and (f) 50-60%. Error bars denote statistical errors. The type B systematic 
uncertainties are represented by the open boxes for the ii2{BBC} and W2{ ZDC-SMD}, and by the solid lines for the V2{2}. The 
gray bands and blue boxes represent the type C systematic uncertainties on the V2{BBC} and U2{ZDC-SMD}, respectively. 
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FIG. 19: The ratio of V2 to U2{BBC} as a function of pt for six centrality bins over the range 0-60% in Au-|-Au collisions 
- 200 GeV. Data symbols are the same as in the Fig. 1181 Error bars denote statistical errors. The solid red lines 



at ./sT 



represent the type B systematic errors on the V2{2]. The blue and yellow bands represent type C systematic uncertainties on 
U2{ZDC-SMD} and V2{2}. 



tween particles from jets affect the cumulant results, but have less influence on U2{BBC}, as explained in Ref. [2; 
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FIG. 20: Comparison of charged hadron V2 at 1 < pr < 1.2 GeV/c as a function of A'part for i;2{BBC} (filled triangles), U2{ZDC- 
SMD} (filled circles) and 1^2 {2} (filled squares) in Au+Au at — 200 GeV. The error bars represent statistical errors. 

The open boxes represent type B systematic uncertainties on W2{BBC} and i;2{ZDC-SMD}. Type B systematic uncertainties 
on 112 {2} are represented by solid red lines. The gray and blue bands represent type C systematic errors on «2{BBC} and 
i;2{ZDC-SMD}, respectively. i'2{2} values are shifted in the x-axis to improve the plot. 
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FIG. 21: Anisotropy parameter V2 as a function of pseudorapidity within the PHENIX central arms using event planes from 
the BBC (filled circles), ZDC-SMD (filled squares), and from the two-particle cumulant method (open triangles) for centrality 
20-40%. The results are shown for three pr bins, which are from top to bottom: 2.0-3.0, 1.2-1.4 and 0.6-0.8 GeV/c. Only 
statistical errors are shown. 



where it was shown that the smaller the rapidity gap be- 
tween the leading particle from a jet and the event plane, 
the greater the V2 of the leading particle of the jet. 

In order to illustrate more clearly the differences be- 
tween the different methods. Fig. shows the ratio of 
W2{ZDC-SMD} and V2{2} to wziBBC}. The results from 
the three methods are comparable in magnitude within 
systematic errors, except for the central and peripheral 
bins where the largest deviations occur. In addition, 
V2{2} and W2{ZDC-SMD} show different behaviors at pT 
> 3 GeV/c, with V2{2} being larger, and U2{ZDC-SMD}, 



smaller than V2{BBC}. 



B. Centrality dependence of V2 

Figure [201 shows the A'part dependence of V2 from dif- 
ferent methods for charged hadrons in the range 1.0 < pt 
< 1.2 GeV/c. V2 is observed to increase with decreasing 
A^part and then decrease slightly for iVpart ^ 75. Note 
that V2 values obtained with the different methods agree 
well within systematic errors for all centralities. This is 
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FIG. 22: (a) Comparison of the U2{ZDC-SMD} from the S-N (filled circles) and ZDC-BBC-CNT subevents (filled squares) as a 
function of pT in the 20-60% centrality range, (b) The same comparison as (a) for the U2{BBC}, where filled triangles and open 
circles represent the v-2 from the S-N and ZDC-BBC-CNT subevents, respectively, (c) Comparison of vi between BBC (filled 
triangles) and ZDC-SMD event planes (filled circles) from the S-N subevent as a function of pr in centrality 20-60%. (d) The 
same comparison as (c) from the ZDC-BBC-CNT subevent, where filled squares and open circles represent the 112 {ZDC-SMD} 
and U2{BBC}, respectively. Error bars denote statistical errors. Open boxes and shaded bands describe the quadratic sum of 
type B and C systematic uncertainties from the S-N and ZDC-BBC-CNT subevents, respectively. 



•pT dependent, as shown in Fig. 1181 



C. Pseudorapidity dependence of V2 

Figure [21] compares the pseudorapidity dependence of 
the W2 of charged hadrons within the r/ range (± 0.35) of 
the PHENIX central arms for different -pT selections. It 
can be observed that vi is constant over the r\ coverage 
of the PHENIX detector and the constancy does not de- 
pend on pt- This is not the case when the V2 is measured 
far from midrapidity where the PHOBOS and STAR col- 
laborations observe a drop in vi for |?7| > 1.0 j47l. |48|. 



VI. DISCUSSION 

A. Effect of CNT event plane resolution 

Figure [55] shows the comparison of -(;2{ZDC-SMD} 
and i'2{BBC} as a function of "pr corrected either by 
the resolution from South-North correlations from the 
same detectors or by the resolution from ZDC-SMD- 
CNT correlations in the 20-60% centrality bin. Fig- 



ures I22f a) and (b) compare the V2 obtained by us- 
ing two different corrections from the South-North and 
ZDC-BBC-CNT subevents for the BBC (a) and ZDC- 
SMD event planes (b). The from the South- 
North subevent is consistent with that from the ZDC- 
BBC-CNT subevent, within systematic uncertainties. 
The small difference between South-North and ZDC- 
BBC-CNT subevents is attributed to the difference be- 
tween the event plane resolution, as shown in Fig. [4] 
Figures [Hfc) and (d) compare W2{ZDC-SMD} with 
W2{BBC} for the South-North (c) and ZDC-BBC-CNT 
subevent (d). The data points in Fig. [5Wc) and (d) are 
the same as in Fig.lSSla) and (b). Figure[22|c) shows that 
V2{ZDC-SMD} is about 10% smaller than wajBBC} for 
the South-North subevent. The ratio of wajZDC-SMD} 
to U2{BBC} is found to be independent of pT except for 
6 < pt < 8 GeV/c. If jets are the dominant source of 
non-flow, one expects its contribution to V2 to become 
larger at higher pT- The constant ratio suggests that the 
non-flow contribution from jets is small and V2 fluctua- 
tions may affect V2{BBC} below px ~ 6 GeV/c since the 
effect of fluctuations is expected to be independent of pr- 
V2{ZDC-SMD} agrees with ?;2{BBC} within systematic 
uncertainties for the ZDC-BBC-CNT subevent as shown 
in Fig. [22yd). The event plane resolution from the ZDC- 
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BBC-CNT subcvents includes the effect of non-flow con- 
tributions and V2 fluctuations since the CNT and BBC 
event planes are sensitive to both effects, though non- 
flow effects especially from jets could be negligible in the 
BBC event plane, as discussed earlier. The consistency 
between V2 from the ZDC-SMD and BBC event planes 
may suggest that W2{ZDC-SMD} becomes sensitive to V2 
fluctuations by the inclusion of the BBC and CNT event 
planes to estimate the resolution. 



B. Comparison with other experiments 
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FIG. 23: (a) Comparison of charged hadron i;2{2} between 
PHENIX (filled squares) and STAR experiments (open stars) 
as a function of px in centrality 20-60%. Solid lines repre- 
sent the quadratic sum of type B and C systematic errors on 
the PHENIX V2{2}. (b) Comparison of charged hadron V2 
from four-particle cumulant D2{4} at STAR (open stars) with 
the PHENIX i;2{BBC} (filled triangles) and waiZDC-SMD} 
(filled circles) as a function of pr in centrality 20-60%. Open 
boxes and shaded bands represent the quadratic sum of type 
B and C systematic errors on the ii2{BBC} and ii2 {ZDC- 
SMD}, respectively. STAR results are taken from Ref. [ioj . 
Systematic errors on the STAR V2 are not plotted, see text 
for more details. 

It is instructive to compare measurements made by dif- 
ferent experiments at RHIC. Figure [53] shows a compari- 
son of the pt dependence of charged hadron V2 in the 20- 
60% centrality range between PHENIX and STAR exper- 
iments [13 . The relative systematic errors on the STAR 
f2{2} and W2{4} measurements range up to 10% for px 
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FIG. 24: Comparison of the PHENIX W2{BBC} (filled tri- 
angles) and U2 {ZDC-SMD} (filled circles) with the STAR V2 
from modified event plane method (open stars) for charged 
hadrons [ij as a function of pr in centrality (a) 10-20%, (b) 
20-30%, and (c) 30-40%. Open boxes and shaded bands rep- 
resent the quadratic sum of type B and C systematic errors 
on U2{BBC} and «2 {ZDC-SMD}, respectively. 



< 1 GeV/c, with the lowest pr bin having the largest 
error ~ 10%, while they are of the order of 1% above 1 
GeV/c m. The z;2{2} from PHENIX is lower than that 
from STAR, but they are comparable within systematic 
uncertainties, as shown in Fig. [231 a). Figure [23jb) com- 
pares U2{BBC} and U2{ZDC-SMD} with W2{4}, obtained 
from four particle cumulants, as measured in STAR. For 
Pt > 2 GcV/c, the STAR W2{4} is systematically smaller 
than the PHENIX event plane V2, while W2{ZDC-SMD} 
is lower than W2{BBC}. However, the three set of mea- 
surements are consistent within systematic errors. The 
order of V2, W2{BBC} > z;2{ZDC-SMD} > z;2l4| could 
be explained by the effect of flow fluctuations [SQ , [sij if 
other non-flow contributions are small. 

Figure [23] shows the comparison of our charged hadron 
V2 from the BBC and ZDC-SMD event planes to V2 from 
a modified event plane method [13], labelled W2{EP2}, 
from the STAR experiment for three centrality bins in the 
range 10-40%. Particles within IA77I < 0.5 around the 
highest Pt particle were excluded for the determination 
of the modified event plane in order to reduce some of 
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the non-flow effects at liigli pr- We find that f2{BBC} 
agrees well with W2{EP2} over the measured pt range, 
whereas ?;2{ZDC-SMD} is generally slightly smaller than 

W2{EP2} . 



VII. SUMMARY 

In summary we have presented PHENIX elliptic flow 
measurements for unidentified charged hadrons from the 
event plane and the two-particle cumulant methods as a 
function of and centrality at midrapidity {\r]\ < 0.35) 
in Au-|-Au collisions at y/s^ = 200 GeV. The first 
harmonic ZDC-SMD event plane is used to measure V2 
and is compared to V2 from the second harmonic BBC 
event plane in order to understand the possible non-flow 
contributions as well as the effect of V2 fluctuations on 
i;2{BBC}. 

The comparison between V2 from two-particle cumu- 
lant and event plane methods shows that they agree 
within systematic errors. However, non-flow effects from 
jet correlations begin to contribute to the two-particle 
cumulant V2, especially for peripheral collisions and at 
high pt- 

In contrast, non-flow effects on i;2{BBC} are very 
small. The measured W2{BBC} values decrease by about 
3% when the central arm event plane is included in the es- 
timate of the BBC reaction plane resolution. This could 
be due to a partial compensation of the non-flow effects 
on the measured V2, though the results of W2{BBC} with 
and without the CNT event plane resolution are consis- 
tent within systematic errors. The strongest evidence 
that non-flow effects on ?J2{BBC} are small comes from 
the observation that i;2{ZDC-SMD} is comparable with 
W2{BBC} within systematic uncertainties in the 0-40% 
centrality range, and are only ^ 5-10% smaller than 
'y2{BBC} for the 40-60% centrality bin. The magnitude 
of this difference could indicate the level at which non- 
flow effects such as jets or the ridge could impact the mea- 
sured flow. However, PHOBOS has observed the ridge 
to be strongest in central collisions [13 where we observe 
that W2{ ZDC-SMD} is comparable with W2{BBC}. For 
collisions that are more peripheral than 40% centrality, 
PHOBOS observes no ridge[30|, so it is unlikely that our 
observation that ?;2{ZDC-SMD} is ~ 5-10% smaller than 
?;2{BBC} for the 40-60% centrality bin is caused by the 
ridge. Moreover, the difference between W2{ ZDC-SMD} 
and ?;2{BBC} is independent of in the measured cen- 
trality range. 

Due to the large pseudorapidity gap between the event 
plane and the particles detected in the central arms 
spectrometer, and the first harmonic event plane from 
directed fiow by spectator neutrons, t;2{ZDC-SMD} is 
considered to provide an unbiased measure of the el- 
liptic fiow. Within systematic uncertainties the mea- 
sured W2{ZDC-SMD} from PHENIX is consistent with 
V2 from the four particle cumulant method measured by 
the STAR experiment in the 20-60% centrality bin, and 



is also consistent with the STAR V2 from a modified event 
plane method in 10-40% centrality bins. These compar- 
isons (1) further demonstrate the validity of the U2{ZDC- 
SMD}, because both STAR results aim to minimize the 
non-flow effects, (2) reinforce the robustness of the BBC 
event plane method at RHIC, and (3) conflrm previous 
studies of the influence of jets on the measured V2 for 
different rapidity gaps. Hence, W2{BBC} can be used to 
infer constraints on the hydrodynamic behavior of heavy- 
ion collisions at RHIC. 
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APPENDIX: DATA TABLES OF V2 

Tables ITvT - IXl show numerical data in the same units as 
plotted in the figures: pt (GeV/c), U2, type A statisti- 
cal error (Tstat, type B systematic error cr^gt and type C 
systematic error uS^st. 
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TABLE IV: V2{2} as a function of pt in centralities 0-10%, 10-20%, 20-30%, 30-40%, 40-50%, and 50-60%. 



Centrality 



(GeV/c) 




V2 




Cstat 




B 
'-^ syst 




(J 
•^syst 


V2{} 


PT (GeV/c) 




V2 




Cstat 




B 
^syst 




Q 





247 





00859 





00014 





00001 





00000 







250 





00898 





00021 





00001 





00000 





347 





01406 





00019 





00004 





00000 







349 





04323 





00030 





00026 





00000 





450 





01882 





00023 





00007 





00000 







448 





06214 





00036 





00053 





00000 





547 





02140 





00027 





00009 





00000 




Q 


548 





07193 





00042 





00071 





00000 





649 





02395 





00031 





00011 





00000 




Q 


648 





08243 





00048 





00093 





00000 





748 





02718 





00036 





00014 





00000 







748 





09401 





00055 





00121 





00000 





847 





03087 





00041 





00018 





00000 




Q 


848 





10533 





00063 





00153 





00000 





949 





03605 





00047 





00024 





00000 







948 





11678 





00071 





00187 





00000 


1 


090 





03950 





00041 





00029 





00000 




1 


092 





12972 





00063 





00231 





00000 


1 


291 





04734 





00053 





00042 





00000 


30-40% 


1 


291 





15059 





00081 





00312 





00000 


1 


490 





05633 





00070 





00059 





00000 


V2{2} 


1 


489 





16955 





00107 





00395 





00000 


1 


689 





06542 





00095 





00080 





00000 




1 


689 





18422 





00147 





00467 





00000 


1 


890 





07148 





00128 





00096 





00000 




1 


891 





19625 





00198 





00529 





00000 


2 


194 





08352 





00128 





00130 





00000 




2 


197 





21718 





00196 





00648 





00000 


2 


698 





09362 





00249 





00164 





00000 




2 


702 





22835 





00369 





00717 





00000 


3 


329 





08866 





00421 





00147 





00000 




3 


338 





22623 





00556 





00704 





00000 


4 


365 





08997 





01134 





00151 





00000 




4 


360 





19059 





01496 





00499 





00000 


5 


376 





07933 





02365 





00118 





00000 




5 


379 





16931 





03256 





00394 





00000 


6 


695 





08701 





02720 





00142 





00000 




6 


628 





16346 





05010 





00367 





00000 


Q 


248 





01089 





00013 





00002 





00000 







250 





00625 





00032 





00001 





00000 





348 





02714 





00018 





00011 





00000 







349 





04611 





00044 





00028 





00000 





449 





03914 





00023 





00023 





00000 







448 





06387 





00054 





00054 





00000 





547 





04592 





00027 





00032 





00000 




Q 


548 





07455 





00062 





00073 





00000 





649 





05281 





00030 





00042 





00000 







648 





08575 





00072 





00097 





00000 





748 





05977 





00035 





00054 





00000 




Q 


748 





09774 





00082 





00126 





00000 





848 





06637 





00040 





00066 





00000 




Q 


848 





11126 





00094 





00163 





00000 





948 





07459 





00045 





00083 





00000 







948 





11974 





00108 





00189 





00000 


1 


092 





08249 





00040 





00102 





00000 




1 


092 





13745 





00095 





00249 





00000 


1 


291 





09506 





00051 





00136 





00000 


40-50% 


1 


291 





15672 





00123 





00324 





00000 


1 


490 





10997 





00067 





00181 





00000 


V2{2} 


1 


489 





17633 





00166 





00410 





00000 


1 


689 





12394 





00090 





00230 





00000 




1 


689 





19315 





00229 





00492 





00000 


1 


891 





13378 





00121 





00268 





00000 




1 


891 





20965 





00309 





00580 





00000 


2 


196 





14881 





00121 





00332 





00000 




2 


199 





21909 





00304 





00633 





00000 


2 


699 





16781 





00232 





00422 





00000 




2 


701 





23572 





00567 





00733 





00000 


3 


328 





16669 





00382 





00417 





00000 




3 


344 





24331 





00808 





00781 





00000 


4 


357 





13468 





01047 





00272 





00000 




4 


346 





26575 





02124 





00932 





00000 


5 


371 





14951 





02244 





00335 





00000 




5 


414 





24613 





03288 





00799 





00000 


6 


587 





11931 





02641 





00214 





00000 




6 


566 





17786 





05097 





00417 





00000 





249 





01127 





00015 





00002 





00000 







251 





01201 





00052 





00002 





00000 





349 





03713 





00022 





00019 





00000 







349 





03575 





00056 





00016 





00000 





448 





05370 





00028 





00040 





00000 




Q 


448 





05111 





00063 





00033 





00000 





548 





06252 





00032 





00054 





00000 




Q 


548 





06256 





00071 





00050 





00000 





648 





07147 





00036 





00070 





00000 







648 





07591 





00080 





00073 





00000 





748 





08144 





00041 





00091 





00000 




Q 


748 





08903 





00091 





00101 





00000 





848 





09118 





00047 





00114 





00000 







848 





09965 





00103 





00126 





00000 





948 





10071 





00053 





00139 





00000 




Q 


948 





11124 





00118 





00157 





00000 


1 


092 





11227 





00047 





00173 





00000 




1 


091 





12340 





00103 





00193 





00000 


1 


291 





12982 





00060 





00232 





00000 


50-60% 


1 


290 





14241 





00133 





00257 





00000 


1 


489 





14786 





00079 





00301 





00000 


V2{2} 


1 


489 





16236 





00178 





00334 





00000 


1 


689 





16113 





00107 





00357 





00000 




1 


689 





17737 





00248 





00399 





00000 


1 


891 





17515 





00145 





00422 





00000 




1 


890 





19295 





00337 





00472 





00000 


2 


196 





19364 





00143 





00515 





00000 




2 


198 





21282 





00330 





00575 





00000 


2 


699 





20931 





00271 





00602 





00000 




2 


700 





22201 





00623 





00625 





00000 


3 


333 





20299 





00430 





00567 





00000 




3 


348 





21980 





00917 





00613 





00000 


4 


356 





19729 





01175 





00535 





00000 




4 


373 





24935 





02292 





00789 





00000 


5 


383 





18635 





02567 





00477 





00000 




5 


452 





36285 





05515 





01671 





00000 


6 


611 





15079 





04839 





00313 





00000 




6 


734 





40554 





08167 





02087 





00000 
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TABLE V: V2{2} as a function of pr in centrality 20-60%. 



K^clLbL all L V 


Pt [yjQ V / cj 




V2 


Cstat 


^ syst 


'^syst 


Pt ((^ey / c) 


V2 


0"stat 


„B 

^ syst 


^ syst 




U.ZOi 


U 


0077© 


A AAA1 1 
U.UUUii 


A nonni 
U.UUUUi 


A AAAAA 
U.UUUUU 


1.489 


A 1 AQQ A 

U.i4oo4 


A AAA K Q 
U.UUUOO 


A AAOAO 

u.uuzyz 


A AAAAA 
U.UUUUU 




0.349 





03793 


0.00016 


0.00019 


0.00000 


1.689 


0.16226 


0.00080 


0.00347 


0.00000 




0.448 





05476 


0.00020 


0.00040 


0.00000 


1.890 


0.17456 


0.00108 


0.00402 


0.00000 




0.548 





06374 


0.00023 


0.00054 


0.00000 


2.198 


0.19027 


0.00106 


0.00478 


0.00000 


20-60% 


0.648 





07303 


0.00026 


0.00070 


0.00000 


2.700 


0.20415 


0.00201 


0.00550 


0.00000 


V2{2} 


0.748 





08283 


0.00030 


0.00091 


0.00000 


3.348 


0.21363 


0.00304 


0.00602 


0.00000 




0.848 





09301 


0.00034 


0.00114 


0.00000 


4.373 


0.19568 


0.00653 


0.00505 


0.00000 




0.948 





10247 


0.00039 


0.00139 


0.00000 


5.452 


0.23823 


0.01494 


0.00749 


0.00000 




1.091 





11444 


0.00034 


0.00173 


0.00000 


6.734 


0.18915 


0.02297 


0.00472 


0.00000 




1.290 





13201 


0.00044 


0.00230 


0.00000 













TABLE VL '!;2{BBC} and t;2{ZDC-SMD} from S-N and ZDC-BBC-CNT subevents as a function of pT in centrality 20-60%. 



Centrality 













S-N subevents 








ZDC-BBC-CNT subevents 




(GeV/c) 




V2 




""stat 




B 
'-'syst 




•-^syst 




V2 




0"stat 




B 
'-'syst 




'-'"syst 





247 





02569 





00009 





00049 





00001 





02486 





00009 





00045 





00001 





348 





04271 





00009 





00016 





00003 





04133 





00010 





00015 





00003 





448 





05587 





00010 





00014 





00006 





05407 





00012 





00013 





00005 





548 





06846 





00011 





00015 





00009 





06625 





00013 





00014 





00008 





648 





08009 





00013 





00015 





00012 





07751 





00015 





00014 





00011 





748 





09123 





00014 





00016 





00015 





08828 





00017 





00015 





00014 





848 





10124 





00016 





00019 





00019 





09798 





00019 





00018 





00018 





948 





11159 





00018 





00017 





00023 





10799 





00021 





00016 





00021 


1 


092 





12439 





00016 





00018 





00029 





12038 





00020 





00017 





00027 


1 


292 





14170 





00020 





00019 





00037 





13713 





00025 





00018 





00035 


1 


492 





15770 





00027 





00027 





00046 





15261 





00031 





00025 





00043 


1 


692 





17244 





00037 





00027 





00055 





16688 





00040 





00026 





00051 


1 


892 





18481 





00050 





00030 





00063 





17885 





00052 





00028 





00059 


2 


200 





19684 





00049 





00029 





00071 





19049 





00052 





00027 





00067 


2 


703 





20803 





00092 





00025 





00080 





20132 





00092 





00023 





00075 


3 


343 





20569 





00141 





00039 





00078 





19905 





00138 





00037 





00073 


4 


381 





17942 





00371 





00066 





00059 





17363 





00360 





00062 





00056 


5 


410 





14862 





00877 





00098 





00041 





14382 





00849 





00092 





00038 


6 


852 





16262 





01770 





00328 





00049 





15738 





01713 





00308 





00046 





247 





02532 





00025 





00047 





00004 





02661 





00035 





00052 





00002 





348 





04002 





00029 





00014 





00010 





04188 





00037 





00015 





00004 





448 





05165 





00032 





00012 





00017 





05395 





00041 





00013 





00007 





548 





06296 





00036 





00013 





00025 





06567 





00046 





00014 





00010 





648 





07433 





00041 





00013 





00035 





07746 





00051 





00014 





00014 





748 





08377 





00046 





00013 





00044 





08730 





00057 





00015 





00017 





848 





09429 





00052 





00017 





00056 





09827 





00065 





00018 





00022 





948 





10365 





00059 





00015 





00067 





10808 





00074 





00016 





00027 


1 


092 





11617 





00053 





00016 





00085 





12065 





00063 





00017 





00033 


1 


292 





13006 





00066 





00016 





00106 





13535 





00081 





00018 





00042 


1 


492 





14367 





00086 





00023 





00129 





14994 





00109 





00024 





00052 


1 


692 





15763 





00115 





00023 





00156 





16504 





00150 





00025 





00062 


1 


892 





17281 





00151 





00026 





00187 





18136 





00203 





00029 





00075 


2 


200 





18031 





00149 





00024 





00204 





18912 





00200 





00026 





00082 


2 


703 





18983 





00263 





00021 





00226 





19998 





00375 





00023 





00092 


3 


343 





18147 





00393 





00030 





00206 





19147 





00576 





00034 





00084 


4 


381 





16102 





01018 





00053 





00162 





17005 





01517 





00059 





00066 


5 


410 





14043 





02402 





00088 





00124 





14833 





03585 





00098 





00050 


6 


852 





12310 





04849 





00188 





00095 





13003 





07240 





00210 





00039 



20-60% 
«2{BBC} 



20-60% 
«2{ZDC-SMD} 
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TABLE VII: ■uaiBBC} from S-N and ZDC-BBC-CNT subevents as a function of in centrality 0-10%, 10-20%, and 20-30%. 



Centrality 






S-N subevents 




ZDC-BBC-CNT subevents 


«2{} 


PT (GeV/c) 




















0.247 


0.01025 


0.00012 


0.00008 


0.00001 


0.00966 


0.00016 


0.00007 


0.00000 




0.348 


0.01868 


0.00014 


0.00003 


0.00002 


0.01762 


0.00025 


0.00003 


0.00002 




0.448 


0.02300 


0.00016 


0.00005 


0.00003 


0.02169 


0.00030 


0.00005 


0.00002 




0.548 


0.02741 


0.00018 


0.00007 


0.00004 


0.02586 


0.00035 


0.00006 


0.00003 




0.648 


0.03174 


0.00020 


0.00007 


0.00005 


0.02993 


0.00041 


0.00006 


0.00005 




0.748 


0.03570 


0.00023 


0.00007 


0.00006 


0.03367 


0.00046 


0.00006 


0.00006 




0.848 


0.03990 


0.00026 


0.00007 


0.00008 


0.03763 


0.00051 


0.00006 


0.00007 




0.948 


0.04428 


0.00029 


0.00008 


0.00010 


0.04176 


0.00057 


0.00007 


0.00009 




1.092 


0.04941 


0.00025 


0.00008 


0.00012 


0.04660 


0.00061 


0.00007 


0.00011 


0-10% 


1.292 


0.05631 


0.00032 


0.00008 


0.00016 


0.05310 


0.00070 


0.00007 


0.00014 


W2{BBC} 


1.492 


0.06349 


0.00042 


0.00008 


0.00020 


0.05988 


0.00082 


0.00007 


0.00018 




1.692 


0.07065 


0.00058 


0.00012 


0.00025 


0.06663 


0.00096 


0.00010 


0.00022 




1.892 


0.07859 


0.00078 


0.00011 


0.00031 


0.07412 


0.00115 


0.00010 


0.00028 




2.200 


0.08557 


0.00078 


0.00009 


0.00037 


0.08070 


0.00121 


0.00008 


0.00033 




2.703 


0.09598 


0.00151 


0.00015 


0.00046 


0.09052 


0.00179 


0.00014 


0.00041 




3.343 


0.09806 


0.00245 


0.00031 


0.00049 


0.09249 


0.00257 


0.00028 


0.00043 




4.381 


0.08795 


0.00699 


0.00089 


0.00039 


0.08295 


0.00667 


0.00079 


0.00035 




5.410 


— 






— 












6.852 


— 


_ 


_ 


— 





_ 





_ 




0.247 


0.01804 


0.00010 


0.00008 


0.00000 


0.01754 


0.00011 


0.00007 


0.00000 




0.348 


0.03095 


0.00011 


0.00008 


0.00001 


0.03008 


0.00015 


0.00008 


0.00001 




0.448 


0.03927 


0.00012 


0.00012 


0.00002 


0.03816 


0.00018 


0.00011 


0.00002 




0.548 


0.04714 


0.00014 


0.00018 


0.00003 


0.04582 


0.00020 


0.00017 


0.00003 




0.648 


0.05480 


0.00015 


0.00016 


0.00004 


0.05326 


0.00023 


0.00015 


0.00004 




0.748 


0.06236 


0.00017 


0.00016 


0.00006 


0.06060 


0.00026 


0.00015 


0.00005 




0.848 


0.06895 


0.00019 


0.00016 


0.00007 


0.06701 


0.00029 


0.00015 


0.00006 




0.948 


0.07647 


0.00022 


0.00017 


0.00008 


0.07432 


0.00033 


0.00016 


0.00008 




1.092 


0.08498 


0.00019 


0.00018 


0.00010 


0.08259 


0.00033 


0.00017 


0.00010 


10-20% 


1.292 


0.09731 


0.00024 


0.00018 


0.00014 


0.09457 


0.00040 


0.00017 


0.00013 


t;2{BBC} 


1.492 


0.10883 


0.00032 


0.00022 


0.00017 


0.10576 


0.00047 


0.00021 


0.00016 




1.692 


0.12204 


0.00044 


0.00021 


0.00021 


0.11860 


0.00058 


0.00020 


0.00020 




1.892 


0.13129 


0.00059 


0.00029 


0.00025 


0.12760 


0.00072 


0.00027 


0.00023 




2.200 


0.14375 


0.00058 


0.00021 


0.00030 


0.13970 


0.00074 


0.00020 


0.00028 




2.703 


0.15569 


0.00112 


0.00023 


0.00035 


0.15130 


0.00120 


0.00022 


0.00033 




3.343 


0.15885 


0.00177 


0.00033 


0.00037 


0.15437 


0.00180 


0.00031 


0.00034 




4.381 


0.13970 


0.00491 


0.00056 


0.00028 


0.13577 


0.00480 


0.00053 


0.00027 




5.410 


0.12763 


0.01194 


0.00101 


0.00024 


0.12403 


0.01161 


0.00095 


0.00022 




6.852 


0.10820 


0.02401 


0.00193 


0.00017 


0.10515 


0.02334 


0.00183 


0.00016 




0.247 


0.02367 


0.00011 


0.00032 


0.00001 


0.02303 


0.00012 


0.00030 


0.00001 




0.348 


0.03981 


0.00012 


0.00014 


0.00002 


0.03874 


0.00015 


0.00014 


0.00002 




0.448 


0.05138 


0.00014 


0.00016 


0.00004 


0.04999 


0.00017 


0.00015 


0.00004 




0.548 


0.06250 


0.00015 


0.00019 


0.00006 


0.06081 


0.00020 


0.00018 


0.00005 




0.648 


0.07276 


0.00017 


0.00020 


0.00008 


0.07080 


0.00023 


0.00019 


0.00007 




0.748 


0.08298 


0.00019 


0.00018 


0.00010 


0.08075 


0.00026 


0.00017 


0.00010 




0.848 


0.09184 


0.00022 


0.00020 


0.00012 


0.08937 


0.00029 


0.00019 


0.00012 




0.948 


0.10139 


0.00024 


0.00020 


0.00015 


0.09866 


0.00032 


0.00019 


0.00014 




1.092 


0.11279 


0.00021 


0.00022 


0.00019 


0.10976 


0.00032 


0.00021 


0.00018 


20-30% 


1.292 


0.12862 


0.00027 


0.00023 


0.00024 


0.12516 


0.00038 


0.00022 


0.00023 


W2{BBC} 


1.492 


0.14459 


0.00036 


0.00029 


0.00031 


0.14070 


0.00046 


0.00027 


0.00029 




1.692 


0.15864 


0.00049 


0.00030 


0.00037 


0.15437 


0.00058 


0.00029 


0.00035 




1.892 


0.17169 


0.00066 


0.00032 


0.00043 


0.16707 


0.00074 


0.00030 


0.00041 




2.200 


0.18437 


0.00065 


0.00032 


0.00050 


0.17941 


0.00075 


0.00030 


0.00047 




2.703 


0.19554 


0.00123 


0.00042 


0.00056 


0.19028 


0.00127 


0.00039 


0.00053 




3.343 


0.19585 


0.00192 


0.00048 


0.00056 


0.19058 


0.00192 


0.00046 


0.00053 




4.381 


0.18189 


0.00521 


0.00088 


0.00049 


0.17700 


0.00509 


0.00083 


0.00046 




5.410 


0.14502 


0.01244 


0.00138 


0.00031 


0.14112 


0.01211 


0.00131 


0.00029 




6.852 


0.15856 


0.02490 


0.00286 


0.00037 


0.15430 


0.02423 


0.00271 


0.00035 
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TABLE VIII: i;2{BBC} from S-N and ZDC-BBC-CNT subevents as a function of pr in centrality 30-40%, 40-50%, and 50-60%. 



Centrality 






S-N subevents 




ZDC-BBC-CNT subevents 


«2{} 


PT (GeV/c) 










V2 










0.247 


0.02733 


0.00015 


0.00064 


0.00001 


0.02643 


0.00016 


0.00059 


0.00001 




0.348 


0.04523 


0.00016 


0.00017 


0.00004 


0.04375 


0.00018 


0.00016 


0.00003 




0.448 


0.05935 


0.00018 


0.00017 


0.00006 


0.05740 


0.00021 


0.00016 


0.00006 




0.548 


0.07263 


0.00020 


0.00016 


0.00010 


0.07024 


0.00024 


0.00015 


0.00009 




0.648 


0.08502 


0.00023 


0.00015 


0.00013 


0.08223 


0.00028 


0.00014 


0.00012 




0.748 


0.09651 


0.00025 


0.00020 


0.00017 


0.09334 


0.00031 


0.00018 


0.00016 




0.848 


0.10742 


0.00029 


0.00019 


0.00021 


0.10390 


0.00035 


0.00018 


0.00020 




0.948 


0.11793 


0.00033 


0.00018 


0.00025 


0.11406 


0.00039 


0.00017 


0.00024 




1.092 


0.13156 


0.00028 


0.00022 


0.00032 


0.12724 


0.00038 


0.00020 


0.00030 


30-40% 


1.292 


0.15004 


0.00036 


0.00019 


0.00041 


0.14512 


0.00046 


0.00018 


0.00038 


z;2{BBC} 


1.492 


0.16604 


0.00048 


0.00030 


0.00050 


0.16059 


0.00057 


0.00028 


0.00047 




1.692 


0.18107 


0.00066 


0.00029 


0.00060 


0.17513 


0.00073 


0.00027 


0.00056 




1.892 


0.19290 


0.00089 


0.00034 


0.00068 


0.18657 


0.00094 


0.00032 


0.00063 




2.200 


0.20640 


0.00088 


0.00035 


0.00078 


0.19962 


0.00094 


0.00032 


0.00073 




2.703 


0.21859 


0.00164 


0.00042 


0.00087 


0.21142 


0.00164 


0.00040 


0.00081 




3.343 


0.21843 


0.00252 


0.00037 


0.00087 


0.21127 


0.00247 


0.00034 


0.00081 




4.381 


0.18342 


0.00662 


0.00101 


0.00061 


0.17740 


0.00641 


0.00095 


0.00057 




5.410 


0.15970 


0.01568 


0.00197 


0.00046 


0.15446 


0.01517 


0.00184 


0.00043 




6.852 


0.18703 


0.03171 


0.00640 


0.00064 


0.18090 


0.03067 


0.00599 


0.00060 




0.247 


0.02840 


0.00024 


0.00071 


0.00002 


0.02735 


0.00024 


0.00066 


0.00002 




0.348 


0.04699 


0.00025 


0.00018 


0.00005 


0.04524 


0.00027 


0.00017 


0.00005 




0.448 


0.06236 


0.00028 


0.00015 


0.00009 


0.06005 


0.00031 


0.00014 


0.00009 




0.548 


0.07757 


0.00031 


0.00015 


0.00014 


0.07469 


0.00035 


0.00014 


0.00013 




0.648 


0.09141 


0.00035 


0.00015 


0.00020 


0.08802 


0.00040 


0.00014 


0.00018 




0.748 


0.10354 


0.00039 


0.00016 


0.00025 


0.09969 


0.00045 


0.00015 


0.00024 




0.848 


0.11530 


0.00044 


0.00019 


0.00032 


0.11102 


0.00051 


0.00017 


0.00029 




0.948 


0.12668 


0.00050 


0.00016 


0.00038 


0.12198 


0.00057 


0.00015 


0.00035 




1.092 


0.14106 


0.00044 


0.00015 


0.00047 


0.13583 


0.00054 


0.00014 


0.00044 


40-50% 


1.292 


0.15967 


0.00056 


0.00019 


0.00061 


0.15374 


0.00066 


0.00017 


0.00056 


t;2{BBC} 


1.492 


0.17584 


0.00075 


0.00025 


0.00074 


0.16932 


0.00083 


0.00023 


0.00068 




1.692 


0.19082 


0.00104 


0.00031 


0.00087 


0.18373 


0.00110 


0.00029 


0.00080 




1.892 


0.20216 


0.00141 


0.00031 


0.00097 


0.19466 


0.00144 


0.00029 


0.00090 




2.200 


0.21274 


0.00138 


0.00031 


0.00108 


0.20485 


0.00142 


0.00029 


0.00100 




2.703 


0.22348 


0.00256 


0.00039 


0.00119 


0.21518 


0.00252 


0.00036 


0.00110 




3.343 


0.22044 


0.00387 


0.00067 


0.00116 


0.21226 


0.00376 


0.00063 


0.00107 




4.381 


0.18665 


0.00994 


0.00094 


0.00083 


0.17973 


0.00958 


0.00087 


0.00077 




5.410 


0.16716 


0.02325 


0.00178 


0.00067 


0.16095 


0.02239 


0.00165 


0.00062 




6.852 


0.15951 


0.04732 


0.00616 


0.00060 


0.15359 


0.04556 


0.00571 


0.00056 




0.247 


0.02767 


0.00043 


0.00056 


0.00003 


0.02604 


0.00042 


0.00050 


0.00003 




0.348 


0.04569 


0.00046 


0.00019 


0.00008 


0.04300 


0.00046 


0.00017 


0.00007 




0.448 


0.06193 


0.00050 


0.00018 


0.00014 


0.05828 


0.00052 


0.00016 


0.00013 




0.548 


0.07654 


0.00056 


0.00014 


0.00022 


0.07203 


0.00060 


0.00013 


0.00019 




0.648 


0.08963 


0.00064 


0.00013 


0.00030 


0.08435 


0.00068 


0.00012 


0.00027 




0.748 


0.10358 


0.00072 


0.00014 


0.00040 


0.09747 


0.00077 


0.00012 


0.00036 




0.848 


0.11362 


0.00082 


0.00020 


0.00048 


0.10692 


0.00087 


0.00018 


0.00043 




0.948 


0.12637 


0.00093 


0.00011 


0.00060 


0.11892 


0.00099 


0.00010 


0.00053 




1.092 


0.14117 


0.00082 


0.00014 


0.00075 


0.13284 


0.00091 


0.00012 


0.00066 


50-60% 


1.292 


0.15953 


0.00105 


0.00020 


0.00095 


0.15013 


0.00114 


0.00017 


0.00085 


W2{BBC} 


1.492 


0.17233 


0.00141 


0.00028 


0.00111 


0.16217 


0.00146 


0.00024 


0.00099 




1.692 


0.18714 


0.00196 


0.00029 


0.00131 


0.17611 


0.00196 


0.00026 


0.00116 




1.892 


0.19757 


0.00266 


0.00054 


0.00146 


0.18592 


0.00260 


0.00047 


0.00130 




2.200 


0.20146 


0.00260 


0.00054 


0.00152 


0.18959 


0.00255 


0.00048 


0.00135 




2.703 


0.21521 


0.00480 


0.00066 


0.00174 


0.20252 


0.00458 


0.00059 


0.00154 




3.343 


0.19757 


0.00712 


0.00083 


0.00146 


0.18593 


0.00674 


0.00074 


0.00130 




4.381 


0.16368 


0.01791 


0.00363 


0.00100 


0.15403 


0.01686 


0.00321 


0.00089 




5.410 


0.11745 


0.04124 


0.00292 


0.00052 


0.11053 


0.03881 


0.00259 


0.00046 




6.852 
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TABLE IX: ualZDC-SMD} from S-N and ZDC-BBC-CNT subevents as a function of pr in centralities 0-10%, 10-20%, and 
20-30%. 



Centrcility 






S-N subevents 




ZDC-BBC-CNT subevents 


V2{} 


Pt (^ije V /cj 




Cstat 


B 
^syst 


c 

^ syst 


V2 


0"stat 


B 
^syst 


c 

''"syst 




0.247 


U.Ulo4z 


0.00114 


0.00013 


0.00005 


0.01723 


0.00158 


0.00022 


0.00004 




0.348 


0.01488 


0.00131 


0.00002 


0.00007 


0.01929 


0.00183 


0.00003 


0.00005 




0.448 




0.00133 


0.00002 


f\ riAAAC 

0.00005 


0.01688 


0.00205 


0.00003 


0.00004 




0.548 


0.02085 


0.00170 


0.00004 


A AAA1 

0.00013 


0.02643 


0.00230 


0.00006 


0.00010 




0.648 


0.01557 


0.00166 


0.00002 


A AAAAT 

0.00007 


0.02132 


0.00256 


0.00003 


0.00006 




n ^ AO 
0. do 


0.02236 


0.00203 


0.00003 


A AAA1 K 

O.OOOio 


0.02928 


0.00289 


0.00005 


0.00012 




0.848 


U.UiooD 


0.00233 


0.00003 


A AAA01 

0.00021 


0.03444 


0.00326 


0.00005 


0.00017 




0.948 


U.U0UI4 


0.00265 


0.00004 


A AAAOT 
0.0002 i 


0.03909 


0.00371 


0.00006 


0.00021 




1.092 


0.04275 


0.00268 


0.00006 


A AAA CT C 

0.00055 


0.04922 


0.00319 


0.00008 


0.00034 


L)-1U70 


1.292 


(J.U382D 


O.UO0O4 


0.UUUU4 


A AAA A A 

0.00044 


A A/l CA1 

0. 04801 


A AA/1 riK 

0.00400 


A AAAAfi 
0.00000 


A AAAQO 
O.OOO02 




1 /IQ9 


u.uoooy 


0.00367 


0.00003 




A A 1 1 ,1 

0.05124 


A A A r /I 

0.00534 


A AAAAI" 

0.00005 


A AAAOT 

0.00037 




1 tiOO 

i.oyz 


n H/i /I no 


r\ r\r\ A 

0.00476 


0.00005 


n AAAftI 

O.OOOoi 


0.06137 


0.00730 


0.00009 


0.00053 




i.oaz 


U.UOoio 


0.00654 


0.00007 


A AA1 OA 
0.00120 


A AO 1" 

0.08583 


A AAAAO 

0.00992 


A AAA1 A 

0.00014 


A AA 1 AO 

0.00103 




z.zUU 


U-UoylU 


0.00672 


0.00006 


A AA1 A Q 

0.00i4o 


A AAOOO 

0.09233 


A AAA A 

0.00989 


A AAA1 1 

0.00011 


A A A 1 1 A 

0.00119 




Z. / Uo 




r\ r\'i 1 

0.01123 


0.00010 


A AA1 QO 

0.00io2 


0.11270 


0.01925 


0.00021 


0.00178 




Q Q y1 Q 

o.Mo 


U.U i 4oi 


0.01667 


0.00018 


A AA1 i^Q 
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TABLE X: U2{ZDC-SMD} from S-N and ZDC-BBC-CNT subevents as a function of in centralities 30-40%, 40-50%, and 
50-60%. 
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